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Abstract 

STSP seeks a pair of pickup and delivery tours in two distinct networks, where the 
two tours are related by LIFO contraints. We address here the problem approxima- 
bility. We notably establish that asymmetric MaxSTSP and MinSTSP12 are APX, 
and propose a heuristic that yields to a 1/2, 3/4 and 3/2 standard approximation 
for respectively Max2STSP, Max2STSP12 and Min2STSP12. 
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1 Introduction 



The multiple Stack Traveling Salesman Probem, STSP, has been recently in- 
troduced in [1] and is very well motivated by the initial request of a trans- 
portation company. Consider a company that owns a vehicle fleet and a depot 
in m distinct cities Ci, . . . , C m . For any pair (Cj, Cji) of cities, the company 
handles orders from Cj suppliers to Cy customers, proceeding as follows: in 
Cj, a single vehicle picks up all the orders that involve a Cj> customer; during 
this tour, the commodities are packed within a single container that consists 
of a given number k of rows; the container is then sent to Cj/, where a local 
vehicle delivers the commodities. The operator thus is faced with two levels of 
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transportation: the local one inside the cities, the long-haul one between the 
cities. STSP models the local level problem; namely, in the context we have 
just described, every pair of cities defines an instance of STSP. The problem 
specificity relies on the way the containers are managed: their rows behave 
like stacks, and no repacking is allowed. The rows of the containers thus are 
subject to LIFO "Last In First Out" constraints. Therefore, the delivery tour 
in Cji must handle at first orders that have been handled at last in Cj. 

Formally, an instance I of STSP considers a number k of stacks (the 
container rows), together with two distinct networks I A = (G A ,d A ), I B = 
(G B ,al B ), where V(G a ) = [n+l] and al a : E(G a ) -> N, a e {A, B}. In both 
networks, vertex represents the depot, whereas vertex % e [n] represents 
the ith order (the ith supplier in I A , the iih customer in I B ). A solution 
consists of a pickup tour T A on J" 4 , a delivery tour T B on I B , and a packing 
V = {P 1 , ■ ■ ■ ,P h } of the commodities into the stacks, where a stack P 13 is 
described as a sequence (ii, . . . , i p p) of commodities (index 1 refers to the bot- 
tom of P 13 ). We denote by < a , a G {A, B} the complete order that T a induces 

— * 

on [n\. i < a i' iff T a handles order i earlier than it handles i' . Similary, we 
denote by < v the partial order that V induces on [n]: i < v %' iff V stacks % 
at a lower position than %' in a same stack P". The LIFO constraints may 
thus be expressed in terms of the three orders < A , < B , <T . Namely, a triple 
(V,f A ,f B ) is consistent iff it satisfies: 

(1) Vi ^ i! e [n], % < v t' (i < A i') A (i > B i') 

The value of a solution (V, T A , T B ) is given by the sum d A (T A ) + d B (T B ) 
of the distance of its two tours. The ultimate goal naturally consists in finding 
the feasible triple (V \T A ,T B ) of optimum value. 

In what follows, the stacks are assumed to have unlimited capacities. The 
most often, input graphs are complete. For TSP {the Traveling Salesman 
Problem), as well as for STSP, A(S)TSP, A(S)TSP and (S)TSPa6 respec- 
tively refer the asymmetric, the metric and the bivaluated cases. Finally, 
/cSTSP considers that k is a universal constant. Since one manipulates both 
maximization and minimization goals, we use notations y, y, opt, opt in- 
stead of >, >, max, min (resp., <, <, min, max) if the goal is to maximize 
(resp., to minimize). Finally, because one jointly considers pickup and delivery 
tours, we respectively denote by E^ 1 = \ (i,i r ) € E}, G~ l = {y,E~ r ), 

d- 1 : E- 1 ->■ N : d-\i,i') = d(i',i) and I~ l = (G -1 , d" 1 ) the reverse of an 
edge set E, a graph G — (V,E), a distance d : E — >■ N, a network / = (G, d). 
Hence, T is a tour of value d(T) on I iffT^ 1 is a tour of value d(T) on J -1 . 



2 Preliminaries 



The problem trivially is NP — hard, from TSP. Indeed, for extremal val- 
ues of the number of stacks, STSP strictly is equivalent to TSP: if STSP 
assumes a single stack, then it reduces to TSP, considering the distance 
d : i — y d A (i,i') + d B (i',i); if on the opposite STSP assumes n stacks, 
then it reduces to two independent TSP (packings that store one commodity 
per stack do not induce any LIFO constraint). The litterature on computa- 
tional aspects of STSP is rather thin: see [1,2] for metaheuristic approachs, 
[3] for an exact approach (that solves 2STSP via the computation of the k 
best tours in I A ,I B ), [4] for an analysis of STSP complexity. We here make 
use of natural connections to TSP and observations made in [4] in order to 
locate STSP within the approximation hierarchy. 

Besides its practical impacts, STSP presents a rather intricate combina- 
torial structure. One the one hand, one could wonder on what part of the 
problem does structure / impact the most the solution / its complexity. It 
appeared in [4] that both the packing and the tour subproblems are tractable: 

Theorem 2.1 ([4]) In arbitrary input graphs: 

(i) Deciding, given a pair (T A ,T B ), whether it admits or not a consistent 
packing V (and to design this packing if the answer is YES) is in P. 

(ii) When A; is a universal constant: given a packing V , computing a pair 
of tours that are optimal for V is in P . 

On the other hand, STSP has strong connections to coloring and schedul- 
ing. Indeed, subproblem (i) reduces to graph coloring in C(T A ,T B ) = {i ^ 
i' £ [n] | i < A i' i < B i'}, whereas subproblem (ii) reduces to scheduling 
under < v (or > v ) precedence constraints. 

3 A first classification of STSP 

Property 1 (i) There exists a standard approximation preserving reduction 
from ASTSP to ATSP that maps a p-approximate solution of ATSP to 
a p/2, p/2 + (a/26), p/2 + (b/ 2a) -approximate solution of ASTSP for 
Max(S)TSP, Max(S)TSPab and Min(S)TSPab, respectively. 

(ii) Case of arbitrary graphs. Any approximate factor p for any restriction 
of STSP also holds for the same restricton of TSP. 

Proof, i. Associate with the instance I = (I A , I B ) of ASTSP the pair (I A , I B ) 



(iii) 

Fig. 1. The three configurations that contradict consistency. 

of instances of ATSP. Compute T a ,a £ {A, B} an approximate tour on I a . 
Pick the tour T A or T B that performs the best among d A (T A ) + (d B )~ 1 (T A ) 1 
(d A )-\f B ) + d A (f B ). ii. Associate with the instance I = (G,d) of TSP the 
instance /' = (7, I' 1 ) of STSP □ 

Proposition 3.1 (Immediate from Property 1) 

(i) Max(A)STSP, MaxA(A)STSP, Max(A)STSP01, 
Min(A)STSP12 are standard- approximable within the following factor: 

MaxSTSP MaxASTSP MaxASTSP MaxAASTSP 

61/162-o(l) ; [8] 1/3, [9] 7/16 - o(l), [7] 31/80, [10] 

MaxSTSPOl MaxASTSPOl MmSTSP12 MmASTSP12 
3/7, [12] 3/8, [11] 11/7, [12] 13/8, [11] 

(ii) MaxASTSPe MaxSNP -hard, [5]; MaxSTSP12e APX-hard, [6]; Min- 
STSP is not 2 p ( n+1 ) standard- approximable for any polynomial p, orP = NP 



(0 



(ii) 



4 A matching based heuristic for the symmetric 2STSP 

In what follows, we say that an edge set A is consistent with a packing V iff 
there exists a completion A' of A into a tour T s.t. T is feasible with respect 
to V. We first provide a characterization of consistency for 2STSP. 

Proposition 4.1 (Due to lack of space, the proof is omitted.) Given a 
packing V , J 13 = {0, ... ,^ + 1} denotes the index set on P@ ' , (3 £ {1, 2}. Arti- 
ficial indexes 0,p^ + l are introduced in order to represent the depot vertex, i.e., 
Zq = i % +1 = 0. Let N be a collection of pairwise disjoint elementary chains on 
[n+l], then A is consistent with V iff it satifies the three following properties 
(configurations that contradict consistency are illustrated in Figure 1): 

(i) No jump: V/3 £ {1,2}, Vj < f £ P, ({if, if,} £ a) 

V (3h<h'eJ 3 ^, {if,i^,... ) i 3 ^,if,}CA) =► j'=j + l 

(ii) No crossing edges: Vj ^ f £ J 1 , Wh ^ h! £ J 2 , {{i),i 2 h }, {i],, i%}} C A 

=> (j <f&h< h!) 

(iii) No way back: Vj £ J\ VA £ J 2 , £ N 
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" ■ c 3 c 4 c 3 c 2 °6 c 5 L 4 

Fig. 2. Possible completions of the approximate packing V: case n even. 

c 2 c 3 c 4 c 5 c 6 c 2 c 3 c 4 c 5 c 6 c 2 c 3 c 4 

^ ' 6 * ~~* ^jB 

c 9 C 8 C 7 c 7 c 8 C 9 C 9 C 8 C 7 C 6 C 5 

j, j' ^ 1, j = j' + 1[2] j, j' j'[2] j' = 1 

Fig. 3. The approximate packing when n is even and p = 1. 

The algorithm relies on a pair (M A , M B ) of optimum weight matchings 
on I A ,I B . Given such a pair, together with some additional edge e*, it aims 
at computing a packing V s.t. N a ,a G {A,B}, is consistent with V, where 
N a coincides with M a (resp., M a U {e*}) when n is odd (resp., even). Let H 
denote the multigraph ([n+1], M A UM ), and C 1 , . . . , C p denote its connected 
components. We assume w.l.o.g., that the depot vertex (i.e., node 0) coincides 
with vertex c\ in C l . If n is odd (iff n+1 is even), then every component C h is 
an elementary cycle {c^, . . . , c£ h , } of even length q h > 2. Otherwise, C h = 
{c^, . . . , c h qh , c^}\{{c^, c^ +1 }} is an elementary chain of even length q h — 1 > 
for a lonely index h G [p]. We define the additional edge e* as follows: 

(2) Vp, e* = arg opt { opt {cf(/)} | / G -F(p)}, where: 

a=A,B 

(3) F(l) = { ({0} U {c\, c\}) x ({0} U {c) +1 , c l J) } \{{0, 0}} 

\h<h')= 



(4) F(p) = U p (h<hl)=l V(C h ) x K(C h ' 



The heuristic is described in Algorithm 1: it first computes M a , a G {A, B} 
(Step (i)) and e* (case n even, Step (ii)); it then builds the approximate 
packing V by considering the components the one after each other, starting 
with C 1 (Step (iv)); it finally computes the best tours T a (V),a G {A, B}, 
with respect to the V (Step (v)). Note that before it operates the packing, 
the algorithm may have in Step (iii) to perform some reindexation of the C h , 
depending on e* (notably, if p > 2, then the two connected components that 
are incident to e* must be stacked consecutively). The even case is illustrated 
in Figures 2 (p > 2) and 3 (p = 1). 

Claim 4.2 (i) Algorithm 1 has polynomial-time complexity. 

(ii) N a ,a G {A,B}, is consistent with the approximate packing V. 

(iii) The solutions returned by the algorithm are 3/2, 3/4 and 1/2 standard ap- 
proximate for Min2STSP12, Max2STSP12 and Max2STSP, respectively; 
these ratios are (asymptotically) tight. 



Algorithm 1. An appproximate algorithm for 2STSP. 

(i) Compute M a an optimum weight matching on I a ,a £ {A,B} 

(ii) Determine the connected components of H in such a way that c\ = 
If n is even Then compute e* = aig opt{opt a=A B {d a (f)} j / G F(p)} 

(iii) If n is even and p > 2 Then consider the components in such a way that: 

(a) If V(e*) n V{C l ) G {0, {0}} Then 3/* € [p]\{l} s.t. e* = {tf qh/2] , c^ 1 } 
Else e* = {cj, c 2 } // /t is tafcen modulo [p], j G {2, . . . , q 1 } depends on c\ 

(b) If q 2 /2Ae* = {cj, c 2 } A 3a, {{0, 4}, {c 2 , c 2 2 }} C M Q Then reverse C 2 

(iv) If n is odd or p > 2: 

j 1 <- (If n= 0[2] A e* = {cj,cf} Then j: [(g 1 - l)/2] + 1 otherwise) 
j 2 <- (If ra= 0[2] A e* = {4,cf} A q 2 = 2 Then g 2 : [g 2 /2] otherwise) 

P 1 *~ ( C 2 1 ■ ■ ■ 1 C jl i C l i • • • > C j2 , cf , . . . , C^ 3 ^ , . . . , , . . . , C^ gP ^, 2 -| ) 

r ^ ^V' • • • ' C jl + l' C g 2' ' ' ' ' S' 2 +l' V' ' ' ' ' IV/21+1 ' • • • ' V' ■ ■ ■ ' C \qP /2]+l> 

Else (thus p = 1, e* = {cj, cj,}, j G {1}U{3, . . . , G {1}U{£+1, n+1}): 
If j,f ^ 1 and j = j' + 1[2] Then P 1 <- (4, . . . , c, 1 ), P 2 <- (4 +lJ . . . , c| +1 ) 
Else If j, j' ^ 1 and j = f[2] Then P 1 <- (c^, . . . , c*), P 2 <- (c] +1 , . . . , c l n+l ) 
Else If j' = 1 Then P 1 <- (c\, . . . , cj), P 2 <- (4 +1 , . . . , cj +1 ) 
Else If j = 1 Then P 1 <- (c^, . . . , cj,^), P 2 <- {c l n+l , cj,) 



(v) Compute T a (V), a G {A, P}, the best feasible tour for V = (P 1 , P 2 ) on P* 

rg' l /2l 1 2 r9 h /2l 12 2 12 

|Y7 2 1 + 1 <? h+1 IW 2 1 + 1 q h+1 q 1 q 2 

Fig. 4. Case n even and p > 2: edge e* induces "no jump" in M Q . 

Proof. Claim (i). Step (v) only requires a 0((n+l) 2 )-time, [4]; for the other 
steps, the fact is either famous, or trivial. Claim (ii). N a ,a G {A, B} triv- 
ially satisfies Properties (ii) and (iii) (resp., (i)) of Propositon 4.1, in any case 
(resp., when n is odd or p = 1). Thus assume n = 0[2] and p > 2. Steps (iii-a) 
and (iv) of the algorithm indicate: 3h G [p],3j G [g 6- ] s.i e* = {c^,cf +1 } and 
e* G P 1 x (P 1 U {0}). Af a possibly links Cj (resp., cf +1 ) to vertices (when 
they exist) cjL 1; cj* +1 (resp., , cP£+i) , with cjL x G P 1 (unless e = {q^c 2 }), 
c h j+l G P 2 , 4 +1 G P 1 and cj#i G P 2 (unless e* = {c\,c\} and q 2 = 2). More- 
over, Steps (iii-b) and (iv) ensure that V satisfies: -.({0, 4} € ^ Q A{c 2 ,c 2 2 } G 
M a ) A (c\,c 2 q2 ) G P 1 x P 2 ). Hence, the chain {cj_ 1; cj, c^ +1 , c^+i} (resp., 



{ c j+i; c j j c i + \ c 2 +1 }) ma y nev er jump c^ +1 (resp., c^^) in P 2 (see Figure 4 
for some illustration). Claim (iii). Let (p*, T A '*, T B -*) (resp., (V,T A ,T B )) 
denote an optimum (resp., the approximate) solution, we denote by OPT 
resp., APX) its value; we observe: 

(i) Va in {A, B}, 3 N' a s.t. T' a = N a U iV' a is a tour, and T' a is consistent 
with V. By optimality of T Q with respect to V, we thus have: 

d a (T a ) h d a (M a ) + d a (N' a ) + l^ [2]d Q (e*), a G {A, P}. 

(ii) If n is odd, then any tour on [n + 1] is the union of two perfect matchings 
on [n + 1]. By optimality of M a , we thus have: 2d a (M a ) y d a (T a >*). 

(iii) If n is even, then V T tour on [n + 1], Ve G T, T\{e} is the union of two 
almost perfect perfect matchings on [n + 1] . Furthemore, by connectivity 
of T a >*, 3e a G T a >* n F(p),a G {A,P}. By optimality of M a , we thus 
have: 2d a (M a ) y d a (T a '*) - d a (e a ). 

From (i), (ii), (iii), we deduce: 

(5) 2APX y OPT + 2 da ( N ' a ) + 1 n=o[2](2 ^ d Q (e*) - ^ d Q (e a )) 

What enables to conclude for Max2STSP, when n is odd. When n is even, 
just observe that e* is optimal on F{p), whereas e a G P(p), a G {A, P}: 

(6) 5Z^(e*) h opt{d a (e*)} h opt{d a (e a )} h l/2^d a {e a ) 

a a 
a a 

For the bivalued case, any edge set S satisfies \S\ < d a (S) < 2\S\, a G 
{A, B}. When n is odd, thus consider \T a, *\ = n + 1 = |iV /a | and conclude: 

(7) 1 /20PT < 2(n + 1) < 2 ^ d a (N' a ) < 4(n + 1) < 20PT 

a 

When n is even, |iV /a | = n. Let Q denote the quantity 2(d A (e*)+d B (e*)) — 
(d A (e A ) + d B (e B )): if the goal is to maximize, assume that Q < 2; then 
d A (e*) + d B (e*) = 2 and d A (e A ) + d B (e B ) > 3, what contradicts the fact that 
e* is optimal. Similary, Q may not reach values 5, 6 when minimizing. 

Tight instances I n , a ,b are defined on K n+ i,n > 1, as: d a {u,v) — a if (v — 
u+2Au = 1[4]) or = u+2Au = 0[4]Aa = A) or = m+2Am = 2[4]Aa = B) 
or (v = u + lAu = l[2]Ao = A) or (d = m + 1A«e 0[2] Act = P), d a (u,v) = b 
otherwise. We consider subfamilies l n ,af> s.t. (a, b) = (1, 0) (resp., (2, 1), (1, 2)) 
for Max2STSP (resp., Max2STSP12, Min2STSP12) and n = 8q- 1 (resp., 8q) 
for the odd (resp., even) case. □ 

Theorem 4.3 (Immediate form Claim 4.2) Mm2STSP12, Max2STSP12 
and Max2STSP respectively are 3/2, 3/4 and 1/2 standard approximable. 



5 Conclusion 



Two main questions arise from the proposed analysis: MinASTP approxima- 
bility, STSP differential approximability. 
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